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I. INTRODUCTION 

The inelastic scattering of photons on (quasi-)free 
charged particles, known also as the Compton scatter¬ 
ing of light, is one of the best studied processes in quan¬ 
tum mechanics. This process demonstrates that light is 
more than a classical wave phenomenon and that quan¬ 
tum theory is required in order to explain the frequency 
shifts as well as the angular and polarization distribu¬ 
tion of the scattered light m- For example, the classi¬ 
cal electro-magnetic theory cannot properly describe the 
frequency shifts at low intensity of the incident photons 
[3], although these shifts are derived quite easily from 
the conservation of the (total) energy and momentum of 
the quantum particles involved in the scattering. In the 
framework of quantum theory, a shift in the wavelength 
of the photons occurs since, for an electron at rest, for 
example, a part of the incident photon energy is trans¬ 
ferred to the recoil of the electron. Therefore, an elastic 
scattering of the light can only be assumed if the photon 
energy is negligible, compared to the rest energy of the 
electrons. This low-energy limit to the Compton scatter¬ 
ing of photons can be described by the non-relativistic 
Schrodinger theory [3]. 

Indeed, derivations of the angular distribution and po¬ 
larization of the Compton scattered light can be found 
in many texts but are made usually for plane-wave pho¬ 
tons and electrons 12 IS]. Such a plane-wave approxima¬ 
tion to the Compton scattering applies if the lateral sizes 
of the incident electron and photon beams are (much) 
larger than their wavelength. In contrast, less attention 
has been paid to the Compton scattering of “twisted” 
beams in which each photon (or electron) carries an non¬ 
zero projection of the orbital angular momentum (OAM) 
along the propagation direction, in addition to the spin 
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angular momentum that is related to the polarization of 
the light EM- Here, we therefore investigate how an 
OAM of the incident radiation affects the angular dis¬ 
tribution and polarization of the (Compton) scattered 
light. In particular, we apply the density matrix for¬ 
malism El to explore the angular distribution and po¬ 
larization of the scattered light and compare this to the 
results of a plane-wave scattering. Calculations of the 
angle-differential Compton cross sections have been per¬ 
formed for the scattering of Bessel beams with different 
opening angles and total angular momenta m. While the 
cross section does not depend on to, it is highly sensi¬ 
tive with regard to the opening angle 0 k of the Bessel 
beams. Moreover, the formulation of the problem within 
the density matrix theory clearly highlights why the re¬ 
sults do not depend on to. This can be explained as 
restriction to those elements of the twisted state pho¬ 
tonic density matrix that are diagonal in the momentum 
quantum numbers due to the spatial symmetries of the 
system. 

This work is structured as follows. In Sec.[TTJ we briefly 
review the non-relativistic description of Compton scat¬ 
tering in the framework of the density matrix theory. 
This includes a short account on the electron-photon in¬ 
teraction in the non-relativistic framework as well as the 
quantization of the photon held. We derive the density 
matrix of the Compton scattered photons in terms of the 
usual plane-wave matrix elements which are weighted by 
the initial photonic density matrix in a plane-wave basis. 
The standard case of plane-wave Compton scattering is 
presented for later reference in Sec. |III| The description 
of twisted photon states is presented in Sec. |IV| Detailed 
calculations have been performed for the angular distri¬ 
bution and polarization of Compton scattered twisted 
photons, and will be explained in comparison with the 
plane-wave case. A short summary is finally given in 
Sec. 0 In the Appendix, in addition, we collect all is¬ 
sues related to the normalization of the plane-wave and 
twisted-wave one-particle states and their density matri¬ 
ces. 
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II. DENSITY MATRIX THEORY FOR 
NON-RELATIVISTIC COMPTON SCATTERING 


We consider the scattering of plane-wave or twisted 
photons on either a beam of free electrons as emitted 
from an electron gun ] 12 | . or on a target material with 
a low work function, much lower than the frequency of 
light, such that the electrons can be considered as quasi- 
free. Twisted light has been produced in a wide range of 
frequencies lamca- If the frequency of the photon w 
in the rest frame of the electron is much smaller than the 
electron rest energy, i.e. if the recoil parameter 


hui 
m e c 2 


« 1 , 


( 1 ) 


we can work in the low-energy limit of non-relativistic 
Compton scattering. We then conveniently work in the 
rest frame of the incident electron, where the theoret¬ 
ical description of the scattering process is much eas¬ 
ier. For an electron beam target the results observed in 
the laboratory frame where the electrons are moving can 
be obtained by just performing a proper Lorentz trans¬ 
formation. For electron beams with low kinetic energy, 
If kin £; 1 keV, the quantitative results in the rest frame 
and the laboratory frame differ very little. The size of the 
target should be larger than the lateral size of the vortex 
light beam PH ESI EH- Throughout the paper we use 
units with h = 1 unless stated otherwise. Moreover, we 
employ Gaussian units, where the fine-structure constant 


A. Density Matrix Formalism 

To describe the angular distribution and polarization 
of the scattered photons it is most convenient to use the 
density-matrix theory pun]]. The density matrix for¬ 
malism has been applied just recently to describe the in¬ 
teraction of twisted light with many-electron atoms and 
ions m ■ In this formalism, the system after the collision 
is described by the final state density operator pj , which 
is related to the density operator of the system in the 
initial state before the scattering, pi, by the scattering 
operator, 

Pf = SpiS ] , (2) 

and where S characterizes the interaction of the particles 
during the collision. 

Before the scattering, the electrons and photons are 
initially independent and uncorrelated. The density op¬ 
erator of the initial state can thus be written as the direct 
product of the electronic (pf) and photonic (pf) opera¬ 
tors [U 


| Pi)(Pi\. The initial photon is described by the initial 
state photonic density operator p 7 = / 5 f h ( 7 ), where 7 
refers to a set of quantum numbers to describe that state. 
Below, we will specify the quantum numbers 7 that are 
needed to represent either a plane-wave or a twisted-wave 
photon. 

Let us now write the final state density operator pf 
in a matrix representation in a plane-wave basis, where 
I/) = \pj,kfAf) abbreviates the final plane-wave elec¬ 
tron and photon states, where p f denotes the final elec¬ 
tron momentum while kf and A f stand for the final 
photon momentum and helicity, respectively. Let us 
also introduce complete sets of initial plane-wave states 
|f) = \p,kX), with f \i)(i\ = 1 (for a detailed discussion of 

i 

the orthonormality and completeness of these plane-wave 
bases see the Appendix), to obtain the density matrix in 
the following form: 


(/|P>l/'> = 

f dkdk' (kX\^\k'X')(p l ,k'A'\S^\f')(f\S\p ll kA) 1 

A,A' 

( 4 ) 


where we employed the fact that the initial electron is in 
a plane wave state as discussed above, and where d k is 
an abbreviation for the properly normalized integration 
measure (for details we refer to the Appendix). This 
equation states that in order to calculate the final state 
density matrix for an arbitrary initial photon state p 7 , 
either plane-wave or a twisted-wave or any other photon 
state, we just need to know the ordinary plane-wave S 
matrix elements to describe the physics. Equation Q 
describes how these plane-wave S matrix elements have 
to be weighted by the elements of the initial photonic 
density matrix in the plane-wave basis (kX\p 7 \k r X'). 

The elements of the S matrix themselves can be rep¬ 
resented in a form 


(f\S\i) = -2m5(Ei-E f )(f\f\i), (5) 


where the delta function ensures the conservation of en¬ 
ergy, i.e. the total energy of the initial state particles Ei 
equals the total energy of the final state particles Ef. 
The matrix elements of the transition operator T can be 
calculated using perturbation theory [3]. 

The normalization of the density matrix is purely con¬ 
ventional and for us it is most convenient to normalize 
it to the total cross section. This is achieved by dividing 
out the flux of incident particles [ 22 ]: 


tr P/ 

TVtr (JPi) ’ 


( 6 ) 


Pi = pf ® pf . (3) 

We describe the initial electron as a plane-wave in a 
pure quantum state | pf) with the density operator pf = 


where T and V are the interaction time and volume, 
respectively. In order to calculate the cross section for 
twisted particles, which are spatially localized, we need a 
convenient definition of flux density operator J in terms 
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of the densities of colliding particles times their relative 
velocity [23]. Since twisted beams are not spatially ho¬ 
mogeneous perpendicular to their propagation direction, 
a proper definition of the cross section according to the 
book of Taylor [33] includes an average over the lateral 
structure of the beam of incident particles. We con¬ 
clude that we should define the cross section by means 
of the spatially averaged density of the initial particles 
(n) := V~ 1 fd 3 xn(x) that is proportional to the trace 
of the initial state density operator ( n) = tr p/V. Thus, 
the operator of the averaged density is just the unity 
operator divided by the quantization volume V (see the 
Appendix), and the flux density operator J in 0 can be 

represented as J = i.e. by the relative velocity, 

the speed of light c, times the operators of the averaged 
particle densities of electrons and photons in the initial 
state. 

Using the basis expansion of the density operator, 
Eq. 0, and by employing the energy conservation in 
0 we find that the cross section 0 contains a factor 
S[E f - Ei)6(E f - E t 0 = 6(E f - Ejfe - E v ). For a 
non-zero contribution to the scattering cross section, the 
total energies of the initial state bases used for the ex¬ 
pansion of the final state density matrix, Eq. 0, need 
to be equal. Finally, we can express the scattering cross 
section as [33] 


a = i i dfcdfc 5{Ei - E f ) 

C LI O'y 

7 / A,A' 

x(k\\h\k'\'){ Pi ,k'\'\f'\f)(f\f\p i ,k\). (7) 

Here we already employed the normalization of the initial 
electron plane-wave states trp® 1 = (p, | p r } = 1. This con¬ 
cludes our general discussion of the density matrix for¬ 
malism and the definition of the cross section. The main 
result of this section is the representation of the final- 
state density matrix 0 of Compton scattered light for 
an incident plane-wave or twisted-wave photon beam— 
characterized by the photonic density operator p 7 —in 
terms of the plane-wave matrix elements for the interac¬ 
tion of electrons and photons. 


B. Interaction Between Photons and Electrons 


Let us now turn our attention on the description of the 
interaction between the electrons and the photons. The 
form of the interaction Hamiltonian, 


~ eAp e 2 A 

-Hint =- h 


m e c 


2 m e c 2 


( 8 ) 


follows from the gauge invariant minimal coupling of the 
electromagnetic field to the free electron Hamiltonian 
[1]. The operator A of the electromagnetic vector po¬ 
tential describes the emission or the absorption of one 
photon [3]. Compton scattering is a two-photon pro¬ 
cess: The incident photon is absorbed by the electron 
while the scattered photon is emitted into some other di¬ 
rection. The one-photon interaction operator A ■ p does 
not contribute to the non-relativistic Compton scattering 
amplitude, because the matrix elements of the electron 
momentum operator p vanish in the rest frame of the 
electron [J- Thus, the Compton scattering amplitude 

can be calculated in first-order perturbation theory by 

~2 

means of the two-photon contribution A to the Hamil¬ 
tonian; the transition matrix elements are just given by 
(/|f|t> = (/|ff int |t). 

The photon field operator A that enters the interac¬ 
tion Hamiltonian, Eq. 0 , can be represented by its mode 
expansion into a circularly polarized plane-wave basis 
Uk\(x) = e zk ' x £kA, where the polarization vector £k\ 
is perpendicular to the wave-vector, fc-e^A = 0. There 
are two independent solutions for each fc, denoted by the 
photon helicity A = ±1. In terms of these plane-wave 
modes the photon field operator is given by 

M x ) = ^dkN k [c kA u kA (x) + cJ. A tifc A («)] , (9) 

A 

where we employ the proper integration measure dfc to 
“count” the basis functions (see the Appendix for de¬ 
tails). The creation operator c^ A creates a normalized 
one-photon plane-wave state from the vacuum |feA) = 
4a|0> that is characterized by its linear momentum 
(wave-vector) k and helicity A. The one-photon states 
are normalized as (feA|fcA) = 1. Moreover, the normal¬ 
ization factor Nk = \j2irc/kV is determined such that 
the energy eigenvalues of the free field Hamiltonian are 
just ojk = c|fc| for one-photon states \kA) gj. It is suf¬ 
ficient to know the above representation of the photon 
field operator in a plane-wave basis because we just need 
to calculate the plane-wave matrix elements by means of 
Eqs. 0 and 0. 


C. The Reduced Density Matrix of Compton 
Scattered Photons 

We are going to investigate the angular distribution 
and polarization of the scattered photons, and we are 
not interested in the electron distribution after the scat¬ 
tering has occurred. We therefore have to calculate a 
reduced density matrix pa'A(&//&/), which depends on 
the direction of the scattered photon kf/kf and its po¬ 
larization state A, by tracing out the unobserved final 
electron states. We obtain the reduced density matrix of 
the Compton scattered photons 
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PA'A(kf/k f ) = 


V 


(27r) 2 ctr p- 


fdkdk Jdp f f dk f k 2 f 5{E f -E i )(k\\p 1 \k'\')M* k , x ,{K')M k x(K), (10) 

A,A' 


in terms of the plane-wave matrix elements A4k\{A) = 
(Pf;kfA\Hi n t\ppkX). The calculation of the plane-wave 
matrix elements for non-relativistic Compton scattering 
can be found in textbooks, e.g. in (4], and we only cite 
here the final result: 

■Mki(A) = - , S( Pi + fc - fc f - p f ) . 


e c 2 L 3 sJUEUTf 


( 11 ) 


Here, S denotes a normalized delta function with the 
property <5(0) = 1. For details we refer to the Appendix. 
Because the reduced density matrix (101 contains the 


product of two plane-wave matrix elements we also get 
two delta functions that ensure the conservation of mo¬ 
mentum. Their product can be reformulated as, 


5{p, + k-pj - k^SiPi + k' - pj - k f ) 

= 5(k' - k)5(p f + kf-pi 


fc), (12) 


where we obtain a factor 5(k' - fc), which consumes one 

of the integrations over the plane-wave bases, dk', in 
Eq. (10). For this reason, only those elements of the ini¬ 


tial state photonic density matrix with k - k! remain in 
Eq. (10). Thus, just the momentum-diagonal elements 


of p 7 contribute to the reduced density matrix of the 
scattered photons and the coherences in the off-diagonal 
elements of the initial photonic density matrix are lost. 
The reason for this behaviour is of course the momen¬ 
tum conservation in the plane-wave matrix elements that 
is related to the spatial homogeneity of the system via 
Noether’s theorem [25]. All particles are described as 
plane-waves, except for the initial photons that are pre¬ 
pared in the hitherto unspecified quantum state p 7 . It is 


the spatial homogeneity of the residual system of incident 
and scattered particles which excludes the interference of 
different momentum components of the initial photonic 
state p-y from the reduced density matrix (10) of the scat¬ 
tered photons. 

The above momentum conservation, together with the 
conservation of the total energy 5(Ef - £)) = 5(ojf + 
p 2 /2m e - oj - p 2 /2m e ), determines the frequency of the 
scattered photons. Recalling that we work in the rest 
frame of the incident electron, p i = 0, we find for the 
frequency of the scattered photons 


w/ = u; 


1 - 


hco 
m e c 5 


(1 -cos6»)+o(- ?l ^] 
V m e c z / 


(13) 


where we temporarily reinstated h. The expression for 
utf in (13) accounts for the well-known frequency red- 
shift of Compton scattering [I], which depends on the 
angle 0 between the initial and the scattered photon mo¬ 
mentum vectors. When starting from a fully relativistic 
QED calculation, we would get the above non-relativistic 
frequency shift in (13) the leading order of an expan¬ 
sion in the small recoil parameter r, Eq. 0- Because of 
Eq. ( p~3| ) , the length of the scattered photon’s wave-vector 
\kf \ = ujf(6)/c is completely determined by its direction. 

In the electric dipole approximation, the momentum of 
the photon and its recoiling effect on the electron is ne¬ 
glected, fc = kf = 0. This is a good approximation when¬ 
ever the recoil parameter, Eq. (jTJ| , is negligibly small. 
The non-relativistic Compton scattering becomes elastic 
within the dipole approximation: u>f = u>. Moreover, the 
dipole approximation coincides with the formal classical 
limit h -»■ 0. Within the dipole approximation we obtain 
as final formula for the elements of the reduced density 
matrix the following expression: 


PA'A(kflkf) 


if 

A,A' 


dfc 


(fcA|p 7 |fcA'} 
tr p 1 


^kfA' ’ ^k\ 


) ( £fc / A ' £fcA ) • 


(14) 


D. Stokes parameters and differential cross section 


In the previous subsection we calculated a suitable ex¬ 
pression for the reduced density matrix of Compton scat¬ 
tered photons, Eq. ( fl4| . We now relate the elements of 
the reduced density matrix to the angular distribution 


of the scattered photons, i.e. to their angular differential 
cross section and the polarization properties. According 
to dUGo], the reduced density matrix can be represented 
by the three Stokes parameters P = (Py, P 2 , P 3 ) via 


PA’A{kf) 


d*l( I + P3 Pi-iP 2 \ 
dn2\Pi + iP2 1 ~ P3 / A / A 


(15) 
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From this representation it is easy to obtain the angular 
differential cross section of Compton scattered photons 
as 


777= E PAA(fc/) = P+1+1 + P-1-1 ■ ( 16 ) 

A=±l 


The Stokes parameters are given by 


r, P+1-1 + P-1+1 

-M. - , 

(17) 

P+1+1 + P-1-1 


„ ip+1-1 ~ ip-1+1 

+ 2 - , 

(18) 

P+1+1+ P-1-1 

„ P+1+1-p-1-1 

+ 3 = 

(19) 


As usual, the Stokes parameters Pi and P 2 represent the 
intensity of light that is linearly polarized under different 
angles with respect to the scattering plane. The scatter¬ 
ing plane is defined by the direction of the incident beam 
of light and by the momentum vector of the scattered 
photon. The Stokes parameter P 3 measures the amount 
of light with circular polarization J20l t2S. Moreover, the 
degree of polarization II is defined as the length of the 
vector P 


n = V p i 2 + P 2 + P i ■ (20) 


where r e = e 2 /(?n e c 2 ) ^ 2.8 fm is the classical electron 
radius. 

To become more specific, we need to specify the scat¬ 
tering geometry, and to express the photon polarization 
vectors in terms of the scattering angles. Let us assume 
that the initial photon beam propagates along the 2 -axis, 
ki/ki = ( 0 , 0 , 1 ) T , while the scattered photon propagates 
into the direction kf/kf = (sin 9 cos p, sin 9 sin tp, cos0) T , 
where the scattering angle 9 and azimuthal angle (p are 
the usual polar and azimuthal angles in spherical coordi¬ 
nates. For the polarization vectors of the scattered pho¬ 
tons we give the explicit representation 


Sk f A = £A { 9 ,if) 


1 

7 ^ 


1 cos 9 cos ip - iA sin ip'' 
cos 9 sin p + iA cos ip , 
^ - sin 9 f 


(23) 


which makes evident that the polarization vector is or- 
thonormalized £* kfA -£k f A' - t> AA / and perpendicular to 
the momentum direction kf ■ £k f A = 0. Because the inci¬ 
dent photon propagates along the 2 -axis, its polarization 
vector is just e kiAi = e Ai (0,0) = (l/N/2,iAi/\/2,0) T . 

The differential Compton cross section as a function 
of the scattering angle is just given by using the repre¬ 
sentations (23) of the photons polarization vectors into 
Eq. (22), and by summing over the final state polariza¬ 
tion according to Eq. (16). This yields 


This concludes our discussion of the density matrix for¬ 
malism. We are now ready to study the angular distri¬ 
bution and the polarization properties of Compton scat¬ 
tered light for both plane-wave and twisted photons. 


III. COMPTON SCATTERING OF 
PLANE-WAVE PHOTONS 


Let us now apply the formalism to the standard case 
of the Compton scattering of plane-wave photons, as 
a starting point for later comparison with the case of 
twisted light. Moreover, this will convince us that we 
normalized the reduced density matrix correctly to ob¬ 
tain the differential cross section by comparing with the 
well known results from the literature. 

We now specify the initial photon state as a plane 
wave with wave-vector k t and in a well defined helic- 
ity state A,;, with the photonic initial density operator 
P 7 = PkiAi = \kiA i )(k i A i \. It has the following represen¬ 
tation in a plane-wave basis: 

{k\\pkiAi \k'X') = S(k - k')S xy S(k - ki)S XAi . (21) 


Using its diagonal elements {k\\pk iAi \kX') 
fi\\'8\Ai,S(k - ki ) into Eq. (14) readily yields the 
reduced density matrix of non-relativistic Compton 
scattering of circularly polarized plane-wave photons in 
the dipole approximation 

PA'A(kf) =r c 2 (e; /A ,-£ feiAl ) (e^A-efcw)* (22) 


^EA = ^(1 + COS 2 0)) 
dU 2 V ’ 


(24) 


where 9 denotes the scattering angle, i.e. the angle be¬ 
tween the wave-vectors of the incident and the scat¬ 


tered light. The result (24) for the angular differential 


cross section of Compton scattered light in the dipole 
approximation is well known and can also be obtained 
by means of classical electrodynamics [3]. An integration 
of Eq. (24) over all directions of the scattered photons 


just yields the well known total Thomson cross section 
a = | 7 rr 2 - 665 mb. This comparison with well-known 
results from the literature shows that we normalized the 
final state density matrix correctly to the cross section. 

Similarly to the cross section, we can obtain explicit 
expressions for the three Stokes parameters of the scat¬ 
tered photons as 


sin 2 9 

1 1 + cos 2 9 

P 2 = 0, 

2A,; cos 9 

P:\ = 


1 + cos 2 9 


(25) 

(26) 
(27) 


As we see from Eqs. (25) - |27| , the scattered photons 
are not necessarily circularly polarized, although the ini¬ 
tial photons were. The ratio of the amount of linearly 
and circularly polarized photons varies with the scatter¬ 
ing angle 9. For instance, under 9 = 90° we have P\ = -1 
and the photons are completely linearly polarized in the 
direction perpendicular to the scattering plane. Never¬ 
theless, the scattered photons are fully polarized, II = 1, 
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independent of the scattering angle. This concludes our 
short review of the analytic results for the differential 
cross section and the Stokes parameters for plane-wave 
photons. We will pick them up in the next section where 
we compare them with the case Compton scattering of 
twisted light. 

IV. COMPTON SCATTERING OF 
TWISTED-WAVE PHOTONS 


where <pk is the azimuthal angle that defines the orien¬ 
tation of one particular vector k(tpk) on the momen¬ 
tum cone. The length of these vectors, k = |fc(</?fc)| = 

+ n 2 , are related to the photon frequency w = ck as 
for plane-waves. 

The normalization factor N tw = \j4n 3 /L z RV that ap¬ 
pears in the definition of the amplitudes ||m A is de¬ 
termined such that the twisted one-particle states are 
orthonormalized in the following way: 


We now turn to the main aspect of our paper: The 
calculation of the angular distribution and the polariza¬ 
tion properties of Compton scattered twisted light. For 
that, we need first to construct the initial photonic den¬ 
sity matrix for twisted photons. 


A. Description of Twisted Photon States and the 
Twisted Photon Density Matrix 


The state of a photon in a Bessel beam that propagates 
along the z-axis, briefly referred to as a Bessel-state or 
twisted photon, is characterized by its longitudinal mo¬ 
mentum K||, i.e. the component of the linear momentum 
along the beam’s propagation axis, the modulus of the 
transverse momentum k ± = |fc ± |, the projection of the 
total angular momentum (TAM) onto the propagation 
axis m, and the photon helicity A 123 ■ A photon in a 
twisted one-particle Bessel-state is, thus, characterized 
by the quantum numbers |q) = |k ± k;||toA). It can be rep¬ 
resented by a coherent superposition of plane-wave states 

EBlMa, 

|Ki«l|mA) = Jdk b KlKfirn (k) |feA), (28) 


where we use the proper integration measure dfc for 
plane-wave states (see the Appendix). From Eq. (281 we 
see that A refers to the helicity of the plane-wave compo¬ 
nents of the twisted photon. The amplitudes b K±K ^ m (k) 
are given by 


^K 1 K^m(^k') — A/)w b{k z ^|| )®Kim(^i) (29) 

and are related to the usually employed transverse am¬ 
plitudes (see e.g. eh da cm) 


2 - 7 T . 



»(*!.) = \ (-*) 


«L 


'<5(|&±| - «±) ■ 


(30) 


Because for a photon in a twisted Bessel state K|| and 
re ± are well-defined, all the momentum vectors k of the su¬ 


perposition (28) are lying on a cone in momentum space 


with fixed opening angle dk = arctan(K ± /ft||). The direc¬ 
tion of the momentum vector k on the cone is undefined, 
and can be parametrized as 



f K ± COS (fiff 

/sin dk cos ipff 

k = k(ifk) = 

k ± sinyjfc 

= k I sin dk sin ifk 


l K ll ) 

\ cos d k j 


(31) 


( k ± «H to A | fiq k|| to'A') 

2t r 2 

= 777 — <5mm'<5AA"5(KL - - k|| ) (32) 

UL Z 

and (k ± K||toA|K ± K|| mA) = 1. This corresponds to the 
normalization to one particle per cylindrical volume V = 
nR 2 L z , where both the radius R and the length L z of 
the cylinder are going to infinity (see the Appendix). 

Let us now construct the density operator for twisted 
photons in the pure quantum state ( f28| ) , together with its 
matrix representation in a plane-wave basis. The latter is 


needed to calculate the reduced density matrix, Eq. (141, 


of Compton scattered twisted light. The normalization 
of the one-photon states implies that the twisted-state 
density operator 


Pi = Pn zH m A = |K ± K||mA)(K ± K||TOA| 


(33) 


has unity trace, tr (p KlK||m A) = 1, i.e. it is normalized to 
an average particle density of “one particle per volume 
V”. The matrix elements of the twisted density operator 
(331 in a plane-wave basis are just given by products of 


the amplitudes b K±K . m a, Eq. (29) 


{k\\p KlKimA \k'X') = <5aa'£>aa b KlK[[m (k) 6« lK||m (fc') 

- S X y6 XA e im ^ k -^ (34) 

and they are diagonal in the helicity quantum numbers. 
We stress that only the momentum-off-diagonal elements 
of the density matrix (34) do depend on the projection 


of total angular momentum m. It enters as the difference 
of the vortex phase factors e %mLpk of the two plane-wave 
components k t k !. On the other hand, the momentum- 
diagonal elements of the above density matrix 


sA| mA \k\') 


= S xy 5 XA ^^6(k z - KiiMClfcJ - k l ) 

V K /_l 


(35) 


which enter the calculation of the reduced density matrix 
of Compton scattered photons, Eq. (141, are completely 


independent of m. Therefore, also the angular distribu¬ 
tion and polarization of the Compton scattered photons 
do not dependent of m. 

A few remarks might be in order why the reduced 
density matrix of the scattered photons is independent 
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from to. In this paper we are interested in the angu¬ 
lar distribution and polarization properties of the Comp¬ 
ton scattered photons. We, therefore, project the final 
state density operator onto a basis of plane-wave states 
that can be observed by a usual detector, which mea¬ 
sures the linear momentum of a photon [29]. Thus, all 
involved particles except for the initial twisted photons 
are described as plane-wave states. From the discus¬ 
sion in the previous section we know that the momen¬ 
tum conservation described by the delta function in the 
plane-wave matrix elements enforces the restriction to 
the momentum-diagonal elements of the initial photonic 
density matrix. Thus, we loose the dependence on to 
because only the momentum-diagonal elements of the 
photonic density matrix (34) that describes the initial 
twisted state contribute due to symmetries of the sys¬ 
tem. 


It is known from previous studies that the scattering 
of twisted particles on spatially homogeneous systems, 
such as plane-waves 1231SQ, impact-parameter averaged 
atomic targets mm, impact-parameter averaged po¬ 
tential scattering [32], leads to angular distributions of 
the scattered particles or fluorescence light that are in¬ 
dependent of to. On the other hand, the coherences of 
the initial state density matrix will play a role for scenar¬ 
ios with a spatial inhomogeneity other than the twisted 
beam. For instance the angular distributions do depend 
on to for the collision of a twisted particle with an in¬ 
homogeneous target, like: a second beam of twisted par¬ 
ticles [331 [33], a localized microscopic target such as a 
single atom [SI [23 EMI, or a quantum dot m- 

Another possibility to recover the coherences in the 
off-diagonal elements of the twisted density matrix is to 
look for the angular momentum of the scattered parti¬ 
cles. In fact, it has been shown in [28l 138] that Compton 
backscattered photons do indeed carry orbital angular 
momentum. In order to access the angular momentum of 
the scattered photons, one needs to determine the final- 
state density operator pj in the the basis of the twisted 
states. This requires a suitable detection operator that 
directly measures the orbital angular momentum of the 
scattered photons 


B. Angular distribution and Stokes parameters for 
the scattering of a twisted photon with well-defined 


TAM 


If we substitute the initial state density matrix (35) 


of the twisted photon into Eq. ([14]), and by performing 
the integration over the plane-wave basis in cylindrical 
coordinates f dk = Idpkdk z dk ± k ± , we obtain the 
reduced density matrix for the Compton scattering of a 
twisted photon (in the dipole approximation Uf = ui) . It 
includes an integration over all plane-wave components 


as described by ipk 
'dip k 


PA'A = r ( 


r d ipk t 

J ~27V 


£ fc/A' ’ ( e fc/A ’ e fc(i£>fc)Ai) i 

(36) 


with k(tpi-) from Eq. (311. This reduced density matrix 
can be directly compared with the corresponding result 
for plane-wave light, Eq. (22). 



0 20 40 60 80 100 120 140 160 180 

Scattering angle 9 [°] 

Figure 1. (Color online) Differential cross section dotw/df! 
for Compton scattering of twisted light as a function of the 
scattering angle 9. Results are shown for various values of the 
momentum cone opening angle 9k of the twisted photons. The 
case 9k = 0 is equivalent to plane-wave photons. Of course, 
the limit 9k -> 90° of large cone opening angles can not be 
reached experimentally. 


If we make use of the explicit representation of the po¬ 


larization vectors (23), the integration over the azimuthal 
angle pk in Eq. (36) can be carried out and yields the dif¬ 


ferential cross section 

= — [(1 + cos 2 9k){ 1 + cos 2 9) + 2 sin 2 9k sin 2 $1 
dO 4 L J 


(37) 


where 9 denotes the scattering angle measured from the 
2-axis, and 9k = arctan« ± /ft|| denotes the opening angle 
of the initial twisted photon beam. As anticipated above, 
the differential cross section is independent of the value 
of the projection of total angular momentum to, but it 
does depend on the momentum cone opening angle 9k- 
Note that in the limit 9k -*■ 0 we recover the plane-wave 
result ( |24| ). 

The results for the differential cross section is shown 
in Fig. [T] as function of the scattering angle 9 for var¬ 
ious cone angles 9k- The red solid curve (9k = 0) cor¬ 
responds to the case of initial plane-wave photons and 
shows the well known symmetric angular distribution 
which is minimal at the scattering angle 9 = 90°, where 
it is just 1/2 of the value at forward or backward scat¬ 
tering. For increasing values of the cone opening angle 
9k > 0 the angular distribution of the scattered photons 
gradually changes and the dip at 9 - 90° becomes less 


















pronounced. For sufficiently large 8k > 8* k the distribu¬ 
tion turns around with the maximum of the angular dis¬ 
tribution occurring at 8 - 90° (e.g. for the blue dotted 
curve in Fig. [l]). This crossover occurs at the “magic 
angle” of 9 k = arccos(l/\/3) « 54.7°, where the angular 
distribution is flat. 

In addition to the angular distribution discussed above, 
the momentum cone opening angle 8k of the twisted pho¬ 
ton beam also influences the polarization properties of 
the Compton scattered light. For incident photons in 
the Bessel state, the Stokes parameters are given by 


Pi = 


(1 - 3cos 2 8k) sin 2 8 


(1 + cos 2 8k ) (1 + cos 2 8) + 2 sin 2 8k sin 2 8 

P 2 = 0, 

4A i cos 8k cos 8 

(1 + cos 2 8k)( 1 + cos 2 8) + 2sin 2 8k sin 2 8 


(38) 

(39) 

(40) 


Similar to the differential cross section, the Stokes pa¬ 
rameters do neither depend on the azimuthal angle of 
the scattered photon, nor on the total angular momen¬ 
tum m of the incident twisted beam, as the whole sce¬ 
nario is cylindrically symmetric. Again the results for 
the plane-wave case are reproduced for 8k = 0. 

The results for the Stokes parameters are shown in 
Fig- m as function of the scattering angle 9 for different 
opening angles 8k■ As discussed above for plane-wave 
photons, the values of Pi and P 3 , which quantify the lin¬ 
ear and circular polarization of the scattered radiation, 
respectively, depend on the scattering angle. Both values 
are sensitive to the cone opening angle 8k of the twisted 
light. For instance, as depicted in Fig. [2] (a), the value of 
Pi at 90° scattering decreases in magnitude for increas¬ 
ing values of 9k, starting from Pi = -1 for the plane-wave 
case (9k = 0). In particular, for the “magic angle” of 
9* k = arccos(l/^) the scattered photons are not linearly 
polarized, since Pi = 0 for all scattering angles. For even 
larger 9k > 9 k the value of Pi is positive which indicates a 
change of the plane of linear polarization of the scattered 
photons which are now (partially) polarized in the scat¬ 
tering plane. Since the sign of P 3 depends on the helicity 
A, = ±1 of the incident photons, we show in Fig. |(b) 
the combination P 3 /Ai instead. The values of P3 grad¬ 
ually decrease for increasing 9k, and approach zero for 
9k -*■ 90°. For twisted light with 9k > 0, the degree of 
polarization is smaller than 1, so the scattered photons 
are not fully polarized anymore. For not too large cone 
opening angles 9k the scattered radiation is depolarized 
the strongest at 9 = 90°. For sufficiently large 9k angular 
dependence of the degree of polarization shows the oppo¬ 
site behavior. In particular for 9k -*■ 90° the scattered ra¬ 
diation becomes completely depolarized for forward and 
backward scattering, while the degree of polarization II 
is nonzero at finite scattering angles, see Fig. [2] (c). 


C. Angular distribution for the scattering of a 
superposition of twisted photons 

As discussed above, the angular distribution of Comp¬ 
ton scattered photons and their polarization does not de¬ 
pend on the value of the total angular momentum m if 
the initial light is prepared in a Bessel state with well- 
defined m; it just depends on the opening angle 9k of the 
beam. 

We now examine the Compton scattering of a coher¬ 
ent superposition of two states with equal longitudinal 
and transverse momenta, k ± and «||, equal helicity A, 
but with two different values of total angular momentum 
m2 > mi. Such a superposition is described by the state 
vector 


I7) = Ci Ik^hTOiA) + c 2 |K 1 K||m 2 A), 


(41) 


where the coefficients fulfil |ci| 2 + |c 2 | 2 = 1 and the state 
|y) is normalized (7I7) = 1. The state vector pT] ) lies on 
the Bloch sphere with the two basis states as the poles 
because |y) is a pure quantum state [Si. The experimental 
generation of such superpositions of twisted states has 
been reported, e.g., in [39] . 

Such superpositions of twisted beams have been con¬ 
sidered previously in theoretical studies of scattering and 
atomic absorption processes, e.g. in Refs. [2H ESI SU]. In 
these previous studies it was seen already that the in¬ 
terference between these two TAM eigenstates leads to 
angular distributions of the scattered particles that de¬ 
pend on the difference Am = m 2 - mi of the total an¬ 
gular momentum values of the two beams. This inter¬ 
ference between the two components of the photon state 
I7}, Eq. (41), can be seen in the off-diagonal elements of 
the density matrix of the initial photon state |y) in the 
plane-wave basis 


(27t) 2 

(A;A| (P'y)m'm ] AjA) — 5 \a tt- &(k± ~ ^1 )^(^j| ~ ^j|) 

tx 1 V 


(« 


c 2 e 


| ^11 ^1 ^2 ^ 

iAm((pfc-7r/2) 


* -iAm(ip k -n/2) ^ 


C 2 ' 


(42) 


The interferences are maximized by choosing coefficients 
c n with equal modulus and a relative phase 6 as, e.g., 
Ci = l/\/2 and c 2 = e lS /\/2. These equal-weighted super¬ 
positions are all lying on the equator of the Bloch sphere, 
and where the relative phase 6 = 0.. ,2tt denotes the lon¬ 
gitude [9j. 

Let us note here that still only the momentum-diagonal 
terms of the photonic density matrix p 7 do contribute to 
the scattering cross section. The superpostion of more 
than one twisted state modifies the distribution of plane- 
wave components on the momentum cone so that they 
are no longer uniformly distributed on the cone as for 
the case of a single Bessel state. Instead, the azimuthal 
distribution of plane wave modes is modulated by the dif¬ 
ference of the total angular momentum of the two beams 
Am. This is described by the off-diagonal elements in the 
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(a) Stokes Parameter Pi (b) Stokes Parameter P 3 (c) Degree of polarization Ei 





Scattering angle 0 [°] Scattering angle 6 [°] Scattering angle 6 [°] 

Figure 2. (Color online) The polarization of Compton scattered photons is characterized by the Stokes parameters Pi (left), 
P 3 (middle) and the degree of polarization II (right) as a function of the scattering angle. Results are shown for beams of 
twisted photons with different momentum cone opening angles 6 k■ As in Fig. [T] the case 6 k = 0 is equivalent to plane-wave 
photons. 


density matrix (42) in the space of the two basis states 


of the superposition (41). 


For such a superposition we readily obtain the differ¬ 
ential cross section in the dipole approximation 


d(7‘ 


2tw 2 ' 

= r 


£ / 77 1 1 + cos ( Am fa* - f) +<*)] 


dn 


x \sA f (0,^)-£Ai (0fc,V5fc)| , (43) 
which can be split into two terms, 


dcr 2 tw d(7 t w dtTint 
+ 


dft 


dfl 




(44) 


The first term, with the “1” in the square brackets, is the 
cross section for the case of a single TAM eigenstate, al¬ 
ready discussed in the previous section, cf. Eq. (37). The 


second term dcri n t/dfl describes the interference between 
the two superimposed TAM eigenstates, and is related to 
the off-diagonal elements of the density matrix in Eq. (42) 
by means of cos (A m(<pk - f) + (5), which describes the 
azimuthal modulation of the density of plane-wave states 
on the momentum cone. 

By performing the integration over the azimuthal angle 
ifk the interference of these eigenstates contributes to the 
cross section 


dcTj, 


dfi 


A- sin (tp + d) sin( 20 fe) sin( 2 d), Am=l, 

-7 cos (2tp + S) sin 2 9k sin 2 9 , Am =2, (45) 


0 , 


Am > 3. 


if Am = 1 or 2 and vanishes otherwise due to dipole selec¬ 
tion rules [21] . When we keep the first order non-dipole 
correction due to the electron recoil in the frequency of 
the scattered photons uif + u>i , Eq. ( fl3| ), we find a non¬ 
vanishing interference contribution also for Am = 3, 


dcr int 

dO 


hujj r;f 
m e c 2 16 


sin(3y> + d) sin 3 9k sin 3 9 , (46) 


which is proportional to the small recoil parameter 
huji/m e c 2 « 1, and is therefore much smaller than the 
interference terms for Am = 1 , 2 . 

Figure [3] displays the angular distributions of Compton 
scattered photons for an incident superposition of two 
twisted Bessel photons for Am = 1,2 and for different 
cone opening angles 9k ■ Obviously, the interference term 
depends on the azimuthal angle <p of the scattered pho¬ 
ton, and not on the scattering angle 9 alone. The disks in 
Fig -0 represent a projection of the forward hemisphere 
of scattered photons, as drawn in the top left panel. The 
superpositions of twisted photons break the axial symme¬ 
try of the initial state. The number of azimuthal modula¬ 
tions of plane-wave states on the momentum cone is just 
Am and the orientation of that pattern is determined 
by the relative phase S , i.e. by the latitude of the state 
I 7 ) on the Bloch sphere. Changing the value of S results 
in a rotation of the distributions of scattered photons in 
Fig -0 with respect to the azimuthal angle p. 

Let us briefly discuss the angular distribution in the 
backward-hemisphere, i.e. for scattering angles 9 > 90°. 
The distributions can be easily read off Eqs. (37) and 
(45). In particular, for the case Am = 1 we obtain the re¬ 
lation dcrjttr -1 )/d^(0, <p) = /dfl(-ir - 6, -p-25). 

For A m = 2, on the other hand, we obtain the symmetry 
relation ip) = dcr!^ 1 ^/dtt(ir - 9, <p). 


V. CONCLUSIONS 


In summary, a theoretical study has been performed 
for the Compton scattered of photons from a Bessel beam 
on electrons in the rest frame of the electrons. In the 
long-wavelength limit of the incident radiation and based 
on the non-relativistic Schrodinger’s equation, the den¬ 
sity matrix theory has been used to analyze the angle- 
differential cross section as function of the scattering an¬ 
gle and for Bessel beams with different cone opening an¬ 
gles and total angular momentum. 

Our formulation of the problem in the density matrix 
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Figure 3. (Color online) Angular differential cross section dcr2tw/df2 of Compton scattered light for initial photons in a 
superposition of twisted waves with Am = 1 (top) and Am = 2 (bottom) for different cone opening angles 9k, and for <5 = 0. 
The disks represent a projection of the forward scattering hemisphere (see the drawing in the top left panel), i.e. the radial 
axis of the disks represents sin 9 for scattering angles 0 < 9 < 90°, and the azimuthal axis represents the azimuthal angle <p of 
the scattered photon. For 9k = 0 (left column), again, the cross section coincides with the results of an incident plane-wave 
photon; they are equal for Am = 1 and Am = 2 and they are independent of the azimuthal angle ip. 


formulation explains why the angular distribution and 
polarization of the scattered photons do not depend on 
the value of the projection of angular momentum: Be¬ 
cause of the symmetry of the system only the those ele¬ 
ments of the incident twisted photon density matrix con¬ 
tribute to the reduced density matrix of the scattered 
photons that are diagonal in the momentum quantum 
numbers. The angular momentum value appears only 
via the difference of the vortex phase factors of the plane- 
wave components of the twisted photon and, therefore, 
vanishes on the diagonal of the density matrix. 

We found completely analytical results for the angu¬ 
lar and polarization distribution of Compton scattered 
photons for the scattering of twisted photons. These dis¬ 
tributions are sensitive with regard to the momentum 
cone opening angle Ok- In particular we observed a de¬ 
polarization of the scattered radiation for large values of 
this cone opening angle. In addition, it was found that 
the angular distributions of the scattered photons for a 
superpositions of twisted photon beams with different m 
differ from the case of a single TAM eigenstate if Am = 1 
or 2 due to dipole selection rules. These differences van¬ 
ish for |A?7i| > 2 in the dipole approximation and remain 
small if non-dipole contributions are taken into account. 

For beams of photons and/or electrons of higher en¬ 
ergy, a relativistic treatment of the electron-photon in¬ 
teraction might be more appropriate instead but should 
not change the central results of this work as long as 
the changes in the recoil parameter are moderate. In 
fact, as long as the recoil parameter r in Eq. (jTJ) is small, 


our results can be directly translated to scenarios of in¬ 
verse Compton scattering by applying a suitable Lorentz 
transformation. For inverse Compton scattering, low- 
frequency photons (e.g. optial laser photons) are scat¬ 
tered on ultra-relativistic electrons, and the backscat- 
tered photons’ frequency is Doppler up-shifted to the x- 
ray regime. 


Appendix: Normalization of quantum states of 
plane-wave and twisted photons 

In this Appendix we discuss all necessary details on 
the normalization of the quantum states that enter the 
calculation of the final-state density matrix in the non- 
relativistic Compton scattering of plane-wave or twisted 
light. The usual normalization in a finite box does not 
work in this case, because the twisted photons are cylin¬ 
drical symmetric modes. Moreover, the states of twisted 
photons are represented as a continuous coherent super¬ 
position of plane waves, while in a finite-sized box the mo¬ 
mentum modes are discrete. Therefore, we need to quan¬ 
tize the modes in an infinite volume V -*■ oo. In order to 
correctly normalize twisted-particle quantum states we 
will explicitly keep all the factors of the formally infinite 
volume. We normalize all one-particle quantum states 
| if), i.e. for the electrons and the plane-wave or twisted 
photons, to unity: = 1. In the following we discuss 

what this implies in detail for the electron and photon 
states, their density operators, spatial wavefunctions, as 
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well as for the spatial probability density. 

1. Plane-wave electron states 

Throughout our paper, all initial and final states of the 
electrons are described as plane waves, i.e. as momentum 
eigenstates p\p) = p\p), where p is the electron momen¬ 
tum operator, and the states are characterized by the 
three quantum numbers of the linear momentum eigen¬ 
value p = (PxiPyiPz)- We require that the one-particle 
states are normalized as 

(p\p) = l, (A.l) 


in the previous subsection. The only small difference 
is that the photon states are characterized also by their 
hclicity A in addition to the linear momentum eigenvalue 
k. We thus orthonormalize the plane-wave one particle 
photon states |feA) as 

(kA\k’A') = 6 AA '5(k-k'). (A. 6 ) 

Moreover, we quantize the photon field, expanding the 
photon field operator 

A(x) = ^dkN k [c kA u kA (x) + cJ.A w feA(®)] , (A.7) 

A 


and, hence, for the orthogonality relation 


(p'Ip) = (^) S(p’-p), (A.2) 


and where the normalization 
with the usual interpretation of lim p /_> p S(p' - 

(L/2tt) 3 mg. 

We may define a regularized delta distribution 


A.l) follows from 


B 
pT = 


Kp'-p) : = (jr) $(p'-p), (A.3) 


into a circularly polarized plane-wave basis u kA {x) = 
e lk ' x £k A , and with the same meaning of dfc as for plane- 
wave electrons. The normalization factor Nk = 'Jiirc/kV 
is chosen in such a way that the free Hamiltonian of the 
photon field is given by 

H 1 = ±-Jd 3 x{E + B) 

= ^dfewfc (cj^CfcA + ^ j , (A. 8 ) 

A 


which has the property lim p ,_> p <5(p' - p) = 1. The prop¬ 
erly normalized integration measure for these plane-wave 
states is given by 

f Sp: 'f(i;) A3p (A ' 4) 

In particular, f dpS(p - p') = 1. These definitions also 
provide the correct way of counting the number of final 
states; their density is just dp. Moreover, the complete¬ 
ness relation for the plane-wave electron states is 

y~dp|p)(p| = 1 - (A.5) 

Because the states |p) are normalized, the trace of the 
density operator p p = |p)(p| is just unity tr(p p ) = 1, 
where we also have to take the integration measure B 
when calculating the trace of the density matrix. 

The position space wave function is just given by 
ipp(x) = (x\p) = e lp ' x /\/V. The particle density is de¬ 
termined as the expectation value of the particle density 
operator n(x) = \x)(x\, as n p (x) = trn(x)p p and just 
yields a constant local particle density of n p {x) = 1/V . 
Therefore one usually says that the plane-waves B 
are normalized to one particle in the infinite quantiza¬ 
tion volume V = L 3 -> oo. 


where E and B are the electric and magnetic field oper¬ 
ators, respectively. 

The one-particle states are generated by the photon 
creation operators from the vacuum (zero-photon) state 

|fcA) = 4 A |0>. (A.9) 


The commutation relations of the photon creation and 
annihilation operators 


[cfcA, cL»,] = 6 AA '6(k - k') 


(A.10) 


completely fix the orthonormalization (A. 6 ). 


The proper measure of the final states now includes 
a sum over the helicity states, thus, the completeness 
relation is given by 


^dk |feA)(fcA| = i . (A. 11 ) 

A 


That means we need to include the sum over the two he- 
licity states in the trace of the photonic density matrix in 
order have tr (pkA) = 1 , and with the same interpretation 
of having one particle per volume V as above. The vector 
potential that corresponds to the plane-wave one-photon 
state can be calculated as 


2. Plane-wave photon states 


A kA (x) := ( 0 |A(a;)|fcA) = N k u k a 


We apply the same normalization to the plane-photon 
states as we employed for the plane-wave electron states 


27TC ik- 

kV kA 


(A.12) 
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3. Twisted-wave photon states 


can also be calculated directly from Eq. (A.15) as 


Twisted photons are cylindrical symmetric modes and 
therefore need to be normalized to a cylindrical volume 
V = itR 2 L z with radius R and height L, along the 2 -axis. 
The twisted photon states that we defined in Eq. (28) are 
orthonormalized in the following sense: 


A' | n 1 n^mA) 

2 - 7 T 2 


RL. 


<5(k[| - K||)$(«4 - K ± )<5 m / m <5A'A , (A. 13) 


which implies that (K 1 Ki|mA|K ± K||?nA) = 1. To prove this 
normalization we make use of the identity for the radial 
delta function 


1 r k j_ 

(n KlK||m A>=^ J dipdzdx 1 x 1 ^j^J^ l (K 1 x 1 ) 

, R 
1 7T r 

=— — J dx 1 x 1 n 1 J 2 l (n 1 x 1 ), (A.17) 

o 

and by approximating the radial integral over x ± for large 
R oo by using the asymptotic expansion of the Bessel 
function for large arguments, J m (x) “ \j2/nx cos(x - 
?7Z7r/2 — 7 t/4 ) |38j[4l], yielding 


n 

/ 


R 

da: ± k ± x x J^j(ac ± Xx) rj — . (A.18) 

7 T 


0 

The vector potential that corresponds to the twisted 
one-photon states is given by 


lim S(k[ - k ± ) = — , (A.14) 

k' ± -*k ± 7 r 


where R is the (infinite) radius of the cylindrical nor¬ 
malization volume. The above identification was proven, 
e.g., in Refs. [22, 35]. Moreover, for the longitudinal 
momentum delta function we employ the usual relation 
<5 (k|| = 0) = L z / 27t, where L z is the height of the quanti¬ 
zation cylinder. 

The probability density that can be attributed 
to the twisted one-particle states n KzK . m k(x) = 
tr (n(x)p KlK . m \) turns out the be not spatially constant, 
but instead is given by 


mA (*^ ) - 2^ p’^rn( K Z x l) j 


(A.15) 


where J m are the Bessel functions of the first kind [-TT] . 
The spatially averaged probability density 


A KlK||mA (aO = (0|A(a;)|«: 1 K||mA) 

/ 27 I" d^Aq < ^ ^ ■ x ± +iK\\z ( \ 

V luL z rJ (2n) 2 kA ' (A 9) 

Except for the different normalization factor in front of 
the integral, this coincides with the vector potential em¬ 
ployed in Ref. m to define the twisted light. It was ar¬ 
gued in [27] that A^^ mA (x) describes a beam with well- 
defined projection of total angular momentum m. The 
product of the amplitude a Kim (fe 1 ) and the plane-wave 
polarization vector £*,a was shown to be an eigenfunc¬ 
tion of the 2 -component of the total angular momentum 
operator J z with the eigenvalue m. 

Performing the momentum integrations in the above 
Fourier integral we obtain for the twisted-wave vector 
potential 


ac iimA 


( x ) = 



x £ (~l) me c mg (A)J m _ mg (K 1 x 1 )e i ^ m ~ m ^ ,fi Tj ms 

m= 0,±1 

(A.20) 


{n) : = y Jd 3 xn(x) = ^trp (A.16) 

which enters the definition of the cross section, Eq. 
is just proportional to the trace of the density opera¬ 
tor because the position eigenstates form a complete ba¬ 
sis. That means, also for twisted photons a normalized 
density operator correspond to one particle per volume 
V = ttL z R 2 . The spatially averaged probability density 


where the sum over m s runs over all possible projec¬ 
tions of the photon spin angular momentum onto the 
^-direction and accounts for the coupling of (the projec¬ 
tions of) orbital angular momentum (me) and spin an¬ 
gular momentum ( m s ) to the total angular momentum 
m = me+m s . This representation of the twisted wave vec¬ 
tor potential employs the unit vectors rj 0 = (0,0,1) T and 
r / ±1 = (1, ±i, 0 ) t /\/2 , and the coefficients Co = -^sin0fc 

and c ± i = |(l±Acos0fe) with the momentum cone open¬ 
ing angle tan 6^ = k ± /k p 
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